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Abstract 

In a model with B — L gauge symmetry, right-handed neutrinos may have exotic local B — L 
charge assignment: two of them with B ~ L = —4 and the other one having B — L = 5. Then, it 
is natural to accommodate the right-handed neutrinos with the same B — L charge in a doublet 
of the discrete 6*3 symmetry, and the third one in a singlet. If the Yukawa interactions involving 
right-handed neutrinos are invariant under 53, the quasi-Dirac neutrino scheme arise naturally in 
this model. However, we will show how in this scheme it is possible to give a value for ^13 in 
agreement with the Daya Bay results. For example the 5*3 symmetry has to be broken in the 
Yukawa interactions involving right-handed charged lepton. 
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I. INTRODUCTION 



Usually it is said that neutrino mass eigenstes may be of the Dirac or Majorana type. 
Notwithstanding, the possibility that the definition of the particle and anti-particle is am- 
biguous was pointed out many years ago from two different motivations. The first one, arises 
when Jauch [1] studying the quantization of spinor fields showed that there are anticom- 
mutation relations which cannot follow from the Schwinger's action principle j2|, and that 
it implies the existence of fermion fields with an intermediate nature between those of the 
Dirac and of the Majorana fields. Jauch showed that a real parameter p appears in the 
anticommutation relation between the spinor field at different space-time points and that it 
may has values in the close interval [0, 1]. The Dirac field corresponds to p = 0, while p = 1 
corresponds to the Majorana field. The cases with < p < 1 correspond to what is now 
called pseudo-Dirac or quasi-Dirac fermion (see below). A Lagrangian formulation of this 
kind of fields was done in Ref. Next, Pauli [3] and others j^, considering the conservation 
of global charges as the lepton number, come to the same conclusion about the existence 
of such a kind of fields. If interactions that violate the conservation of these charges are 
allowed, the neutrino field which enter in the weak interactions are of the form au + bu'^, 
with |ap + |6p = 1. Moreover, Pauli pointed out the theoretical possibility that processes 
like the neutrinoless double beta decay may have cross sections with values between zero 
and the theoretical maximum value. The latter point was rediscovery in the context of a 

n 

gauge model in Ref. p. 

In modern gauge theories, since there are several neutrino fiavors, these sort of fields can 
be realized in several ways. They occur when two Majorana neutrinos are mass degenerated 
and have opposite parity so they are equivalent to one Dirac neutrino [7^. However, the 
degenerescence would be remove by the weak interactions, when quantum corrections are 
taken into account the would be a Dirac fermion will split into two Majorana neutrinos 
with different masses. But, if the mass splitting is small this Dirac neutrino will became a 
pseudo-Dirac when one of the mass generated neutrinos is active and the other is sterile N, M 

rinn 

or pseudo-Dirac if both of them are active neutrinos [g, 00, 111. 

□ 

More recently, a quasi-Dirac neutrino scheme was proposed in Ref. [12] . This model uses 
the 5*3 discrete symmetry to generate, at the tree level, the so called tribimaximal mixing 
matrix in the lepton sector. This was also implemented in a model whith exotic right-handed 
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neutrinos i.e., they carry non-convencional U{1)b-l charges [13|]. The aim of the present 
paper is to illustrate the main characteristics of the scheme and its limitations to produce a 
^13 7^ in agreement with the recent Daya Bay results IJ]. Solar neutrino data may impose 
too strong constraints that it is not possible to generated through quantum corrections the 
observed 6^13. 

The outline of this paper is as follows: In Sec. |TT]we review the model in which the scheme 
is implemented. In Sec. Illll we consider the mechanism for generating charged lepton masses 
and in Sec. |IV]we show how the observed 6^13 could be generated. The last section is dedicated 
to our conclusions and some remarks. 



II. MODELS WITH NON-IDENTICAL RIGHT-HANDED NEUTRINOS 



The possibility that all neutrino flavor eigenstates are part Dirac and part Majorana, 
as a consequence of a S3 symmetry, has been put forward recently 12|]: two mass eigen- 
states are Majorana fields and one is a Dirac field at the tree level. This scheme called 
bimodal/schizophrenic is just an interesting example of a quasi-Dirac neutrino based on the 
5*3 symmetry. It is this symmetry that allow to distinguish among lepton generations (dou- 
blets and right-handed neutrinos). However, if all right-handed neutrinos have the same 
quantum number their separation in S3 irreducible representations is arbitrary. This would 
not be the case if the model has two right-handed neutrinos having a different quantum 
number from that of the third one. It happens in the model of Ref. |l3| because of the 
gauged B — L symmetry. In fact, the quasi-Dirac scheme was already implemented in a 
model ISj with local B — L symmetry which has a different scalar content with respect to 



that of Ref. 



13| 



Here we will be concerned on how the non-zero ^^13 may be obtained in the scheme which 
is also consistent with the Daya Bay results [l^. First, let us briefly review the main feature 
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of the model which has the following gauge symmetries: 



SU{3)c ® SU{2)l ® U{1)y' ® U{1)b-l 
SUi3)c®SUi2)L^U{l)Y 

SU{3)c ® (1) 

where Y' is chosen in order to obtain the hypercharge Y of the standard model, given by 
Y = Y' + {B — L). Here (0) and ($) denote one (or several) singlets and doublets of 
SU{2)l, respectively. Thus, in this case, the charge operator is given by 

^ = h + l[Y' + {B-L)]. (2) 

The anomaly cancelation is also implemented if, for quarks, charged leptons and active 
left-handed neutrinos the local B — L charges are as the usual ones, but for right-handed 
neutrinos this charge is, instead of the usual assigment, B — L = —4 for two of them and 
B — L = 5 for the third one. For this reason we call them exotic right-handed neutrinos. 
Thus, these neutrinos are, because of the local B — L and Y'[= —{B — L)], charges, naturally 
split in the two irreducible representation of 5*3: 3 = 2 © 1. After the breaking of the 
electroweak gauge symmetry as is shown in Eq. (|2]), the usual global U{1)b and U{1) l appear 
as accidental symmetries as in the standard model. Also, because of the non-standard local 
charges, the lepton sector has its own scalar sector and no large hierarchy in the Yukawa 
couplings is necessary. In this paper, as we said before, we will consider more details of the 



model of Ref. [15| and suggest two possibilities of how a non-zero 61^ can be obtained. 

The right-handed neutrinos with B — L = —4 are in a doublet of 6*3, say {rieR, ^t-r, may be 
consider heavy, and the singlet say n^R, is the light right-handed neutrino. It is also assumed. 



as in the model of Refs. 



12 



15 



16| (here we called it case (a)), D = (Li, L3) = (l/-\/6)(2Le — 
L^ — Lt), (1/v^)(L^ — Lt-)) transforms as a doublet and As,= L2 = (l/v^)(Le + L^ + L^) as 



a singlet. The scalar sector differs slightly from that of ISj: two scalar doublets with weak 
hypercharge Y = —1 are denoted by $1^2 = {y^i 2 ^12)'^ singlet of 5*3. The mixing angles 
in the (rigij, UrR) sector have been absorbed in h2 and /13. There are also scalar singlets of 



SU{2). For the quantum numbers of these fields see Ref. [l5|. We denote {^i{v2)) = Vi{v2). 
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Therefore, after integrating the heavy degrees of freedom, the effective Yukawa interac- 
tions that give neutrino masses are 



hl 



hi 



rrir, 



[{Ll)n [LiL^] + ^m)n n] [L^l^] + H.c. (3) 



(4) 



From the Yukawa interactions in Eq. ([3]), the neutrino mass matrix in the basis 

f mi \ 

mo 

mg 

\ TTlD / 

which has the eigenvalues mi, m^), ms, m^). At the tree level, there are four massive Majo- 
rana neutrinos, two of them, an active and the sterile neutrino, are mass degenerated and 
correspond to a (quais)Dirac neutrino. 



On the other hand, the neutrino mass matrix (x')^mX' written in the active neutrino 
basis x' = ^'l + {^'lY where iV£ = {veL, ^fj.L, ^t)J,, is of the form 

\ 

(5) 



/ 2 

3 



mi 



11, ms 1 
3 6 ~'~ 2mi 6 



2mi 



v 



1 1 

3 6 



ma 1 
2mi 6 



^ 2mi / 



The eigenvalues of the mass matrix are (mi,0,m3) and will be denoted by mf 



M 



mi. 



m^^ = 0, and m^^ = m^. The massive Majorana neutrinos are Ui and u^, while ud has no 
Majorana mass at tree level, and we have definded mi = /i|w|/m„^ and rris = /igfl/m^^. 
The matrix in Eq. ([5]) is a consequence of the 5*3 symmetry 17| . This matrix is diagonalized 
by the tribimaximal matrix 18| : 

/ /I ^ ^ 



V = Utb^' 



2 J_ 

3 

1 1 ]_ 

]_ j_ 1 



(6) 



where Q' = diag(e*'', 1, 1). The mass eigenstates basis denoted by A^'^, = {uil, i^dl, ^^sl)^ 
are related to the flavor basis as N'^^ = UtbNl- Notice that as a result we have only one 
Majorana phase because one of the eigenvalues of the matrix ([5]) is zero. 

We will understand better this situation if we consider the 4x4 mass matrix, x^Xi where 
Xi = + {NiiY, and Nn = [vil v^l i^sLn^^iV- The eigenvalues of the matrix in Eq. 



are mi, rriD, 771^,171 £, where mi, ma are as given above, and ttid = hiVi/\/2. The respective 



eigenvectors are 



M 



Xl 



1,0,0,0), z/f = (l/v^)(0,X2,0,X4) 



.M 



X3 = (0,0,l,0) and 



i/-\/2)(0, X2, 0, — X4), i-e., we have two Majorana neutrinos, and a Dirac 



neutrino formed by two mass degenerate Majorana neutrinos vd = 1^2 + i^t ■ The matrix 
in Eq. (jlj) has a conserved lepton number induced by the transformation acting only on 
the fields X2 ^~^^X2 and Xi ~^ e~*^X4) implying the Dirac character of vn, which is not 
conserved by the weak interactions, hence we have a quasi-Dirac neutrino formed by an 
active and a sterile neutrino (according to the notation discussed in Sec. [T]). 

Next, let us also consider, from the Yukawa interaction in Eq. (|3]), a 4 x 4 mass term for 
the neutrinos, x' M^^x' 1 but now in the basis x'i = ^[l + {^[lY where N'^j^ = (u^ 
and M'^^ given by: 



M 



mi 





2 


1 




1 


mr> 


\ 


3 


3 




3 






1 


1 _l_ ms 
6 ~'~ 2mi 


1 


_ ma 


rrio 






3 


6 


2mi 


\/3mi 






1 


i _ "^3 


1 


_|_ -mz 
~'~ 2mi 


mn 






3 


6 2m\ 


6 


\/3mi 




I 


rriD 


ran 




rriB 





/ 








\/3m.l 



(7) 



where mi and are defined as in the matrix in ([5]) and m/) as above. The eigenvalues of 
the matrix ([7]), denoted by mf^'^e*'', mg^, m*^ 



, Tfi^ , are 



m 



M 



M 

nil, ^2 



m 



M 



M 

rriD, mg 



(8) 



and we see that the two Majorana masses m^ and are the same as in the case of the 
matrix in Eq. (jS]). The four Majorana massive neutrinos have been split as follows: two of 
them have different masses and are purely Majorana fermions; the other two, which are mass 
degenerated, fuse to form a Dirac massive neutrino. All of this is at tree level. After the 
breaking of the gauge B — L symmetry nothing protects neutrinos to gain small Majorana 
masses by quantum corrections, in particular, the left- and right- components of the would 
be Dirac neutrino, V£,. However, it may be rather small and this neutrino will continue to 
be, for all practical purposes, a Dirac neutrino (see below). 
The mass matrix in ([7]) is diagonalized by the matrix 

^ 1/v^ -l/V6^ 

-i/Ve i/Ve -1/V2 -i/Ve 
-1/V6 1/V6 i/v^ -1/V6 
1/V2 l/v^ 



\ 



(9) 



and Q = diag(e*^, 1, 1, i), which can be rewritten as the following matrix product 16| 



V\ 



(a) 



2 J_ 

3 V3 

1 1 

1 1 



4, 



1 

1 
V2 



\ 



\ 







/ 1 \ 



V 



75 



1 

V2 



10 
VO 73 71 / 



1]. 



(10) 



The relation of the mass eigenstates Xi = + {Ni^Y, with A^'^ = (z/i z/2 z/3 z/4)i with the 
flavor eigenstates x' is given by x = V\(a)X'^ where ^|(a) is given in Eq. or f ITU]) . The 
respective Majorana mass eigenstates fields are i/f^ 



J^iL + {j^ilY and z/|^ = u^l + {v^l 



which have masses mi = mic^P and ms, respectively, and the other two Majorana fields 



V. 



M 



^2L+{t-'2lY^ ^'^d 1/4 = i^iL+ii^iLY form a Dirac field, ud = z/2 +zi/4 , with a Dirac mass 



term, m^,. We define V2udl = ^21 - i^AL = {niY)l and \/2{v%)r = (2/2)^ + i(z/4; 



n 



Note that the mass eigenstates ui and z/3, when written on the flavor basis, do not have the 
contribution of the fourth neutrino n^/j. This is a prediction of the model at tree level, since 
it is the rotation matrix in ( ITOj) which determines this outcome. 
With the masses in (Ml we obtain 



Am2i 



\Aml^\ 



m 



D 



mi 



/)2 



vl 



"'3^2 



2m? 



6 



\Aml^\ 



m. 



m. 



(11) 



Experimentally, Am^^ = (7.59±0.20) x 10"^ eV^ Am|2 = (2.43±0.13) x lO'^eV^ (l9|. Just 
for an illustration, these values can be fitted by choosing, for the normal hierarchy (here and 
below all parameters with dimension of mass are in eV), hi = 0.01, /i2 = 0.025, /i3 = 0.71, 
Vl = 1.85, f 2 = 10^, and m„^ = m^^ = 10^^, the neutrino masses are ttld = 0.0107, mi = 
0.0063, and m^ = 0.0504. We can also obtain the inverted mass hierarchy, for instance by 
choosing hi = 0.0045, /i2 = 0.07, /is = 0.025, vi = 1.9155, t;2 = 10^ and m„^ = m„^ = 10", 
the neutrino masses are m£, = 0.0505, mi = 0.0498, and m^ = 0.0062. Notice that it is a 
prediction of the model the existence of a light sterile neutrino, n^ij. 

In fact, it is possible to choose different representations from that we have called case (a), 
which leads to the effective interactions in Eq. ([3]), in such a way that the neutrino mass 
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term written as (NMN) becomes 

f mi \ 
ma 
mD 
\ niD J 



(b), or 



f mz)\ 
m2 
mg 
\mD / 



(c). 



(12) 



The model of Q, Q corresponds to the case (a), the cases (b) and (c) arise if we define 
b): in the right-handed neutrinos sector rirR is now the singlet and {neR^n^Rj the doublet, 
and Ls = and D = (Li,L2). c): n^R is the singlet and {n^R.riT-R) the doublet, and in 
Ls = Li and D = (L2, L^). A similar analysis to that doing for the matrix in Eq. (jl]), follows 
for matrices in Eq. ( fT2|) . 

Instead of ( JTOl) we have for the cases (b) and (c), respectively: 



V\ 



2 
3 

1 

Ve 



73O 

^ 



\ 



V6 V2 V3 ^ 

1 



/ 1 \ 
10 

71 



\0 



(13) 



V2 V2 



with Q as in Eq. f lTOj) . and 





1 


[2 

V 3 





o\ 




1 








1 






1 

v/3 


1 


1 










1 












1 


1 


1 













1 







v 











V 


v 


1 








1 

V2 


I 



(14) 



where Q." = (l,e*^l,^). 

As we said before, the tribimaximal matrix diagonalized the neutrino mass matrix at 
leading order. Corrections to that matrix may arise from quantum loop corrections [20] (or 
evolution with the renormalization group equations 21|) and/or the mixing in the charge 
lepton sector. In fact, global neutrino data analysis had already suggested that ^13 7^ |22| . 
and a first evidence that this is the case has been obtained from the observation, at the 2.5a 



level, of the appearance — )■ z/g [23|, |2^. More recently, the Daya Bay results are more 
conclusive, at the 5.2a level tjiey found sin2 2^i3 = 0.0920.016(stot)0.005(s|/st) Q, which 



implies 0.13 < sin ^13 < 0.17 



25| . Finite quantum corrections to the mass matrix in Eq. ([5]) 



may be consider but they are strongly suppressed by solar neutrino data because of the 
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active to sterile neutrino oscillation present in this model [26|, l27|. It is for this reason that 
the would be Dirac neutrino is for practical proposes a Dirac fermion. This implies a Vpmns 
matrix of the tribimaxiaml type. The only way to obtain a realistic form of this matrix is 
to have a nondiagonal mass matrix in the charged lepton sector and using Vpj\/Ar5 = UjU^. 
Here, If^ means the 3x3 submatrices in Eqs. (ITU]) . f ll3l) and f[T^ . Next, we have to consider 
how it is possible, in the context of the present model, to obtain a realistic Vpmns matrix. 
Before doing it, let us consider the mass matrix of the charged lepton sector. 



III. THE CHARGED LEPTON SECTOR 



Let us suppose that the 5*3 symmetry also constrains the Yukawa interactions in the 
charged leptons. In this case, to obtain the charged lepton mass matrix we add three extra 
scalar doublets with Y = +1, denoted ^e,^iJ.,^T, and using the complex representation for 
5*3 determined by the matrix 28 1 

/i 1 i\ 

1 cu cu'^ 
1 u"^ u 



(15) 



which change from the irreducible basis 1 and 2 to the reducible one 3 without changing the 
product rules. In this basis 3l = (L^, L^, Lr), 3r = (cr, fiR, tr) and 85 = (<l>e, $r)- 
The charged lepton sector the Yukawa interactions can be written using the four 5*3 singlets 
formed by the direct products of three triplets {xi,yi,Zi), i = 1,2,3: xiyiZi+X2y2Z2 + X3y3Z3, 
xiy2Z3 + xiy3Z2, X2yiZ3 + X3yiZ2, and X2?/3^i + X3y2Zi: 



-Cl = G[Le<i>eeR + L^^^^iR + Lr<^rTR\ + HLe{fiR<^r + T/J^^) 
+ F{L^<l>r + Lr<l>f,)eR + I{L^Tr + LrfiR)<l>e + B.C. 



(16) 



From (11 6p we obtain the most general mass matrix for charged leptons (here {ip^) = vj \/2): 



Ml 



1 
71 



^ Gv, HVr Hv^ ^ 



FVr GVf, IVe 



(17) 



Using the following values for the parameters (dimensional parameters in MeV): Ve = 9.97, 
= 2038.56, = 34924.8, G = 0.05 and H = 1.95 x 10-^ F = 0.35 x 10'^, I = 1.92 x lO'^, 



we obtain rrie = 0.4894, = 102.155 and rrir = 1746.24. The mass matrices are a prediction 
of the model, hence they are vahd at the energies at which all symmetries of the model are 
realized, i.e., at the electroweak scale. For this reason, we use the running lepton masses 
at yU = Mz which values were taken from Ref. {29]. The f/| unitary matrix diagonalize the 
Hermitian MiMj matrix, where M; is the mass matrix in (fT7|) . We obtain for the values of 
the parameter above: 

^ 0.9976 -0.0698 -0.0002 
-0.0698 -0.9976 0.0007 
0.0002 0.0007 0.9999 

We have consider several set of values for the parameters in the Yukawa interactions ([1 
and all of them give similar values for the entries of the matrix as in (ITS]) . 



U 



( 



\ 



\ 



J 



(18) 



IV. THE PMNS MIXING MATRIX 



Notwithstanding, the full PMNS matrix also includes rotations in the charged lepton 
sector, i.e., it is defined as Vpmns = U^l^tb- Here, Utb denotes any of the 3x3 submatrices 
in Eqs. (HDD, (USD and dS]). 

Here, we compute the leptonic mixing matrix, Vpmns 1 for the three cases in Eqs. ffTOl) . 
(fl3|l and (fT4|) . obtaining. 



PMNS I (a) 



7875 





6144 





0473 


4617 





5375 





7056 


4081 





5775 





7070 



\ 



(19) 



and 



I VpMNS I 



\VpMNs\{c) 



( 



5358 





8429 





0495 


6158 





3500 





7059 


5777 





4087 





7066 



\ 



8429 





0492 





5358 


3500 





7049 





6158 


4087 





7076 





5777 



(20) 



(21) 
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respectively. From data we have 25| : 



PMNS\ exp 



\ 



\ 



(22) 



0.78 - 0.845 0.52 - 0.61 0.13 - 017 
0.40 - 0.58 0.39 - 0.65 0.57- 0.80 
0.19 - 0.43 0.53 - 0.74 0.59 - 0.81 

Comparing the three Vpmns in Eqs. f|T9|) - fl2T]) with the matrix in Eq. fj22|) . we see that there 
is no agreement with the measured Vpmns- 

There are a few ways to turn around this trouble. On the one hand, since the neutrino 
and charged lepton masses are at the Z-peak 29| , it means that the running of these values 



to the low energies have to be done still 2l| in the context of the present model. On the 
other hand, radiactive corrections will induce a non-zero 6i^. Notwithstanding, it implies 
that the corrections to the neutrino mass matrix are not to small, but since in this case 
there are active to sterile neutrino oscillations the mass corrections are strongly constrained 
by solar data {2^ in such a way that no realistic value for ^13 arises. 

Another possible way to overcome this difficulty is by considering the S3 symmetry to 
be break in the charged lepton sector. Thus, a S3 non-invariant Yukawa interactions given 
mass to the charged leptons is 



(23) 



where is an arbitrary 3x3 matrix and H denotes the Higgs boson of the standard model 
which gives mass to the quarks too. If we solve, simultaneously the following equations 



Vpmns = U'^Utb, ^U[Y\YyiUiy = diag(m,^ mj, m^). 



(24) 



we obtain that: in case (a) has solutions with complex Yukawa parameters, case (6) has no 
solution and, case (c) has solution with real Yukawa parameters. Here we show only the 
latter case: with 

^0.9114 0.8346 0.6669 ^ 
0.1356 -0.4644 -0.5304 . (25) 
0.3240 0.5457 0.4440 



Ui\{c) 

which implies, from (^^, that 



8.98 1.84 -21.01 
27.886 18.51 33.97 
-26.68 -37.22 -27.14 



(26) 
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and we obtain 



PMNS\{c) 



^ 0.79 0.56 0.13^ 
0.53 0.65 0.71 
0.33 0.58 0.70 



(27) 



which in agreement with 



V. CONCLUSIONS 



The quasi-Dirac neutrino scheme of Refs. 12|, |l5| is interesting in its own. Although there 



is a hght sterile neutrino, its spectrum is not of the form "3+1", since the fourth massive 
neutrino is quasi degenerated with one of the active neutrinos. Hence, the effect of the extra 
mass square difference might appear only in neutrinos coming from long distance sources 



3Q| 



like supernovas and not as a solution to some possible neutrino anomalies 

The mass matrix in equation f ll7p is the best mass matrix that can be obtained without 
violating the symmetry and contribute to the value of ^13 7^ 0. However, without violating 
the 5*3 symmetry in the charged lepton Yukawa interactions we cannot explain the Daya Bay 
data, since the right-handed charged lepton do not have exotic values for the gauge B — L 
symmetry, the 5*3 symmetry is not naturally incorporated in that interactions. Hence, the 
validity of S3 symmetry only in the neutrino sector can be a prediction of the model that 
made the B — L symmetry a gauge symmetry and exotic right-handed neutrinos. On the 
other hand, if the charged lepton interactions violate the symmetry 6*3, we can get a general 
mass matrix, as that in Eq. ( 123|) and, as we have shown, it is possible to fit a realistic PMNS 
mixing matrix in agreement with the Daya Bay's results. 

Some final remarks are in order. 1) If a Majorana mass for the sterile neutrino, say n^R, is 
allowed the 44 entry in the mass matrix ([7]) is non-zero. Notwithstanding, the Majorana mass 



for the n^R arises only from non-renormalizable interactions since the operator {rin^J'^n^ji is 



not allowed by the B — L attribution of the model and the symmetries, see Ref. 15| for the 



respective quantum numbers. Hence, the operator with the lower dimension generating this 



mass term is one of dimension seven: AM('l'sM^'^20i(2)0a.M^)(n^j:j)'^n^R. In this case M = 
^mVsj,jV2UiVx/ J^^ ■ Just for illustration, if \m ~ 0(1), -ui ~ A, M Vgj^^V2Vxl ^ is rather 
small using v^^^ ~ 100 GeV, V2 = 10"^ GeV and A = 1 TeV we have M ^ 10-^^{vx/GeV). 
Moreover, is not necessarily large since it is not the responsible for breaking the B — L 
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gauge symmetry. We see that the M value is rather small and satisfies the constraint from 
solar data [27|. 

2) The present model involves Higgs scalar doublets which couple mainly to leptons and 
have small VEVs. Moreover, the model has also scalar singlets and some of them may have 
VEVs lower than the TeV scale. Tiny VEVs have been proposed before 3l|] and after 32 1 



the models in Ref. 



13 



15|. Then, the question concerning on the stability of the tree level 



VEVs arises. This has been done recently 



33| in the context of the particular 2DHM of 



Ref. 32|. This question is rather model dependent and a similar analysis in the context of 
the present model will be considered elsewhere. 

3) The model in Ref. 13[ was proposed just as a new solution to the anomaly cancellation 
when i? — L is a local symmetry, and it has interesting features by its own. For instance, 
beside implementing the bimodal scheme without fine tuning in the neutrino Yukawa inter- 
actions, it is a model which also implement naturally the features of the so called lep tophilic 
two Higgs doublet model (L2DHM) 3^ and the neutrino specific 2HDM one 32, l35|. In 
the latter models those Higgs doublets were introduced ad hoc. In general the supersym- 
metric versions of this sort of models has interesting features in accelerator physics 



and in cosmology 37|. Recently, a model with quasi-Dirac and Majorana neutrinos in the 
context of supersymmetric standard model with the extra symmetries 5*4 ® {Z^Y has been 



proposed 



38| . However, we would like to stress that all these features, for instance dou- 



blets given masses just for neutrinos, arise naturally when we consider the anomaly free 
B — L gauge symmetry and they are not assumed ad hoc. This avoids extreme fine tuning 
in the Yukawa coupling of the Dirac neutrino as hi = 0.01(0.0045) in the normal (inverted) 



'12 



hierarchy instead hi ~ 10" as m 

4) Finally, we stress that the existence of such a very light sterile neutrino is a prediction 
of the model, and it is not motivated by possible anomalies observed in neutrino experi- 
ments 30||. Since it is (almost) mass degenerate with one of the active neutrinos its effect 
may be only observed in extragalactic neutrinos 16|. 
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Appendix A: charged leptons 



The charged lepton masses in Ref. [12| is generated by a dimension five operator and it 
is diagonal at tree level. Thus the Vpmns is just the matrix that diagonalized the neutrino 
mass matrix and this is just the tribimaximal one. To obtain those authors introduce three 
gauge singlet scalars cXe, (J^ and these fields and the right-handed charged letons transform 
like the left-handed doublets and for avoiding a general mixing it is necessary to impos Z„ 
symmetries: Zn,e ® Zn^^ ® Zn,T in such a way that the right-handed leptons transfom as o;^^ ^ 
and the singlets scalars as oj^^^t- This case has been rouled out by recent neutrino data. 

Case 1. The right-handed components trasfoorm as es = cr, = {^r,tji) and the 
Yukawa coupling is given by 

-Ci = KUesH + h^riLoeoVH + H.c. (Al) 

where H is the usual SM Higgs doublet. This case is not favored because the ^ and r are 
mass degenerated and the matrix which diagonalize the mass matrix is the tribimaximal, 
hence Vpmns = U^p = 1. 
Case 2. Now 



(A2) 

with the Yukawa interactions as in Eq. f lAip . In this case the matrix which diagonalize the 
charged leptons mass matrix is again the the tribimaximal and = (2/3)mr. 

Case 3. The right-handed components of the charged leptons transform as in Case 1 but 
we introduce three Higgs doublets, $£, $r, transforming under S3 as a singlet Hs = 
and a doublet Hd = ($^, $,-). The Yukawa interction is 

- Ci = KL2esHs + h2[L2eD]2HD]i + h3[[Lr)Hs]2eD]i + /i4[[-^D^^D]2eD]i + h[[Lr)Hs]2eD]i 
+ K[LDHD]ies + H.c. (A3) 

In this case we have that DetM' = 0, hence the electron remains massless at tree level. 

Case 4. Introduce Higgs scalars as in Case 3 and all right-handed charged lepton transform 
as singlet under 5*3. The Yukawa interactions is 

-Ci = {hlsHs + h',[LDHD\i)kR + H.c. (A4) 
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In this case we obtain as in the Case 1, the trimaximal matrix in the charged lepton sector 
and two leptons remains massless at tree leveL 

Case 5. Now the e/j,/i/j,r/j transform as in Eq. (lA2p and the three Higgs doublets as 
in Case 3 and the Yukawa interactions are also given in Eq. flA3p . This case is difficult to 
analyse analitically but numerical calculations indicate altought we can fit the three charged 
lepton masses the mixing matrix is given by 



and the PMNS matrix 



^ -0.382 0.621 0.684^ 
0.360 0.582 0.730 
0.851 0.525 0.001 



^ 0.51 0.06 0.86^ 



(A5) 



v 



\VpMNs\ = 0.06 0.99 0.04 (A6) 
0.86 0.03 0.52 ^ 
does not fit the experimental values, see Eq. (!22|) . 

All these cases allow a diagonal mass matrix in the charged lepton sector if extra Z, 
symmetries are added as in Ref. [l2|. 
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